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ABSTRACT 

The abundance of iron is measured from emission line complexes at 6.65 keV 
(Fe line) and 8 keV (Fe/Ni line) in RHESSI X-ray spectra during solar flares. 
Spectra during long-duration flares with steady declines were selected, with an 
isothermal assumption and improved data analysis methods over previous work. 
Two spectral fitting models give comparable results, viz. an iron abundance that 
is lower than previous coronal values but higher than photospheric values. In 
the preferred method, the estimated Fe abundance is A(Fe) = 7.91 ± 0.10 (on 
a logarithmic scale, with A(H) = 12), or 2.6 ± 0.6 times the photospheric Fe 
abundance. Our estimate is based on a detailed analysis of 1,898 spectra taken 
during 20 flares. No variation from flare to flare is indicated. This argues for a 
fractionation mechanism similar to quiet-Sun plasma. The new value of A(Fe) 
has important implications for radiation loss curves, which are estimated. 

Subject headings: Sun: abundances — Sun: corona — Sun: flares — Sun: X-rays, 
gamma rays — line: identification 
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INTRODUCTION 



The solar abundance of iron remains an important parameter and topic in solar 
physics. Iron is the most abundant of all elements with Z > 14, and is a large contributor 
to the radiation loss at coronal temp eratures. Recent dete rminations of the photospheric 



abun dance give A(Fe) = 7.50 ± 0.04 ( lAsplund et al. 



20091 ) and 7.52 ± 0.06 (ICaffau et al. 



201 ll ) (on a logarithmic scale wh ere A(H) = 12), in n ear-agreement with the meteoritic 



abundance, A(Fe) = 7.45 ± 0.01 (ILodders et al. 



20091 ) . The iron abundance in the corona 



has been determined from X-ray or extreme ultraviolet Fe emission lines, formed by 
collisional excitation. As the excitation rates are a function of electron temperature 
T e , the thermal structure of the emitting coronal plasma must be modeled for correct 
interpretation of line fluxes, and ionization fractions and excitation r ate co efficients must 



be known. Examp l es of t he coro nal Fe abundance in clude 



Fludra fe Schmelz 



(119991 ) (7.65), 



White et al. 



(120001 ) (8.19), and 



Parkinson! (119771) (A(Fe ) 



Dennis fe Starr! (12008! ) 



7 65). 



(7.86), i.e. enhancement factors over the photospheric value of between 1.4 and 4.9. This 
large range of abundance determ inations may indicate time variations in the coronal 



abundance ( iSylwester et al. 



19841 ) or measurement uncertainties of ~ 0.5 in A(Fe). The 



Fe abundance from solar energetic particles (SEPs) in the interpla netary medium varies 



by large factors; a baseline list of abundances for gradual events ( iReames 



19951 ) gives 



A(Fe) = 7.93, or 2.65 times photospheric. Evidence for systema t ic dif ferenc es betwee n 



photo spheric and cor onal abundances has been cited by 



(119921 ). According to 



Meverl ( 119851 ) and 



Feldman 



Feldman fe Laming! (l2000f ) elements with first ionization potential 



FIP < 10 eV like Fe have coronal abundances enhanced by factors of 4, apart from 
low-altitude flares and energetic spray-like events for which the coronal and photospheric 
abundances are equal. Various models have been put forward to explain the fractionation 
and its dependence on FIP; they generally involve a mechanism that separates ions and 
neutral atoms in the chromosphere where low-FIP elements are partly ionized but high-FIP 
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elements are neutral. The mecha nisms include magnetic fields carrying ions rising into the 



corona as active regions develop (IHenoux 



in the acceleration of Alfven waves ( 



Laming 



199 8J), a nd the presence of a ponderomotive force 



2009|) 



Flare spectra in the photon energy range ~ 3 keV to 17 MeV have been obtained from 
the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESST) since its launch on 
2002 February 5, allowing analysis of thermal spectra (range ~ 6-30 keV) with ~ 1 keV 
spectral resolution (FWHM). This ra nge includes th ermal continuum emission (free-free and 



free-bound) and two line complexes ( IPhillips 



20041 ) . The Fe line complex at ~ 6.65 keV is 



made up of Fe XXV lines and dielectronic satellit es of Fe XXIV and lower stages, with min or 



contributions from the Fe XXVI Lyman-a lines (IFeldman et al. 



1980 



Tanaka et al. 



1982 



Lemen et al. 



19841 ). emitted over a broad temperature range (> 10-100 MK). The weaker 
"Fe/Ni line" line complex at ~ 8 keV consists of higher-excitation lines (Is 2 — lsnp, n ^ 3) 
of Fe XXV and Fe XXIV satellites, with a few-percent contribution from Ni XXVI-Ni XXVIII 
lines. Line-to-continuum ratios can be derived from RHESSI spectra, and from these the 
iron abundance estimated using models for the thermal structure (isothermal or simple 
functions for the temperature distribution). This was do ne for RHESSI spec tra taken 



20061 ) using fluxes 



during various phases of 27 flares between 2002 and 2005 (IPhillips et al. 
of the Fe line complex at 6.65 keV expressed as the equivalent width of nearby continuum. 
The observed equivalent widths were found to follow a dependence on T e , derived from 
the energy dependence of the continuum emission, ap proximately equal to the th eoretical 



dependence with an Fe abundance somewhat less than 



Feldman fe Laming ( )2000l )'s coronal 



value, A(Fe) = 8.10. In this work, this analysis is carried further. First, we have chosen 
only spectra during the gradual phases of flares with the RHESSI thin attenuators in place, 
and secondly we have used advances in the analysis software enabling RHESSI spectra 



to be better interpreted, including 



Young et al. 



2003 



Dere et al 



the use of the chianti code (v. 6) (IDere et al. 



1997 



20091 ) with latest atomic data instead of the earlier mekal 
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code (IMewe et al 



19851 ). Also, instrumental effects can be adjusted in the fitting process. 
By these means, we derive an estimate of the Fe abundance with much reduced uncertainties 
compared with previous work. 



2. SPECTRAL ANALYSIS AND OBSERVATIONS 
2.1. The RHESSI Spectrometer 



Smith et al. 



(120021 ) so only an 



The RHESSI spectrometer has been described by 
outline is given here. Nine cryogenically cooled hyper-pure germanium detectors make up 
the spectrometer, each with a 1-cm-thick front segment which senses low-energy X-rays. 
Modulation collimators in front of each detector result in time-modulated counting rates 
as the spacecraft rotates, which are unscrambled with imaging software to form images. 
The X-ray emission is viewed through and partially absorbed by beryllium windows and 
aluminized Mylar insulation layers. To avoid detector saturation and to reduce pulse 
pile-up, sets of aluminum attenuators are moved in front of the detectors at increasing X-ray 
flux levels. Since instrumental effects like pile-up, slight changes in energy calibration at 
high count rates, and variable background rates are different for each detector, we chose to 
analyze spectra from the detector with best energy resolution rather than summing coun ts 



from multiple detectors to improve statistics. Following our earlier work (IPhillips et al. 



20061 ). we selected flare spectra during times of slowly decaying emission in the Al state 
(thin attenuators in place) when the flare plasma is most nearly isothermal. Spectral fits 
were done over a range from ~ 5.7 keV (the count rates at lower energies are dominated 
by K-escape events in the Al and A3 attenuator states) up to 20-40 keV, depending on 
the emission at higher energies relative to the background spectrum. Thick-target X-ray 
continuum from a power-law electron spectrum was included when necessary to achieve an 
acceptable fit to the data at the higher energies. The energy bins for the spectral fits were 
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those | keV-wide bins used by the on-board pulse-height analyzer. 

2.2. Data Selection 

Twenty flares having a steady decline of X-ray emission were observed by RHESSI 
between 2002 and 2005 in its Al attenuator state during the ~ 60-minute solar-viewing 
part of its orbit. For two long-duration flares (on 2002 July 20/21 and July 26/27), the 
decline could be followed for up to four orbits. Table CD gives the dates and GOES classes 
with time intervals over which analysis of at least the Fe-line complex was possible, one 
interval per flare except for the 2002 July 20/21 and July 26/27 flares; numbers indicate 
flares, letters time intervals for the 2002 July 20/21 and 26/27 flares. Also listed are the 
number of spectra iVf u n, and N COQt analyzed in each interval for each of the chianti_full 
and chianti_cont emission models, and details of the estimated Fe abundance which will 
be discussed in Sectional 



Table 1. Time Intervals Analyzed during RHESSI Flares 



Flare Date RHESSI Time GOES CHIANTI_full Analysis CHiANTi_CONT Analysis 

Range (UT) Class Mean value A(Fe) A(Fc) 

iVf u u ofF A(Fe) N cont Fe-Line Fe/Ni-Line 

Complex Complex 



1.. 


2002 


Mar 10 


22: 


:57- 


23 


:50 


M2 


.3 


9 


0.54 ±0.04 


7. 


.83 


± 





.03 


186 


7 


.91 ±0. 


.03 


8.02 ±0.07 


2.. 


2002 


Apr 15 


00: 


:15 - 


00 


:40 


M3 


.7 


59 


0.56 ±0.19 


7. 


.85 


± 





.13 


88 


7. 


.86 ±0. 


.09 


8.02 ±0.10 


3.. 


2002 


May 31 


00: 


:13 - 


00 


:57 


M2 


.4 




indeterminate 6 












149 


7. 


.93 ±0. 


.13 


8.04 ±0.14 


4 .. 


2002 


Jun 1 


03 


:52 - 


04: 


:04 


Ml 


.6 


17 


0.93 ±0.55 


8. 


.07 


± 





.20 


36 


8. 


.01 ±0. 


,09 


indeterminate 


5a.. 


2002 


Jul 20/21 


22 


:29- 


23: 


:28 


X3. 


3 


21 


0.52 ±0.07 


7. 


.81 


± 





.06 


90 


7. 


.80 ±0. 


,08 


7.87 ±0.20 


5b.. 






00 


:06- 


01: 


:03 






205 


0.53 ±0.22 


7. 


.82 


± 





.16 


172 


7. 


.89 ±0. 


,12 


8.07 ±0.16 


5c. 






01 


:42 - 


02: 


:40 








indeterminate 












172 


7. 


.87 ±0. 


,20 


indeterminate 


5d.. 






03 


:24- 


04 


:16 








indeterminate 












86 


8. 


.03 ±0. 


,21 


indeterminate 


6a.. 


2002 


Jul 26/27 


23 


:01 - 


00 


:00 


M4 


.6 


185 


0.67±0.12 


7. 


.93 


± 





.07 


188 


7. 


.90 ±0. 


.07 


7.97 ±0.12 


6b.. 






00 


:37- 


01 


:36 






59 


0.66 ±0.09 


7. 


.92 


± 





.06 


59 


7. 


.99 ±0. 


.09 


8.09 ±0.10 


6c.. 






02 


:14- 


03 


:13 








indeterminate 












118 


8. 


.04 ±0. 


,13 


8.09 ±0.14 


7.. 


2002 


Jul 29 


10 


:50- 


11: 


:26 


M4 


.7 


2 


0.64 ±0.09 


7. 


.91 


± 





.05 


79 


7. 


.89 ±0. 


.05 


7.91 ±0.07 


8.. 


2002 


Oct 4 


05 


:41 - 


05 


:56 


M4 


.0 


33 


indeterminate 












26 


7. 


.91 ±0. 


,05 


8.03 ±0.12 


9.. 


2002 


Dec 2 


19 


:23- 


19 


:32 


C9. 


6 


23 


0.73 ±0.23 


7. 


.96 


± 





.12 


22 


8. 


.00 ±0. 


.03 


8.05 ±0.10 


10.. 


2002 


Dec 17/18 


23 


:35 - 


01 


:01 


Ml 


.6 


71 


0.43 ±0.25 


7. 


.73 


± 





.20 


29 


7. 


.84 ±0. 


.15 


indeterminate 
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2.3. Spectral Analysis 

For each time interval, data and detector response matrix (DRM) files were extracted 
and read by OSPEX (Object Spectral Executive), an object-oriented IDL program recently 
substantially modified. The non-solar background spectrum was determined from the 
night-time parts of the RHESSI orbit. An isothermal fitting function was chosen to 
model the continuum and line emission, with goodness of fit determined by the reduced 
chi-squared x^ ed , calculated as the weighted sum of the squares of the differences between 
background-subtracted count rates in each energy bin and the predicted count rates 
computed by folding the assumed incident photon spectrum through the DRM. The weights 
were determined from the predicted rates assuming Poisson statistics with zero systematic 
uncertainties. A graphical user interface allowed least-squares spectra with normalized 
residuals to be viewed (see Figures [1] and [2]) . One form of the fitting function we chose, 
unavailable previously, is a thermal spectrum (vth) calculated from chianti, including 
all lines and free-free and free-bound continua. The abundances of individual elements, 
most especially Fe, may be varied independently by a factor F from a part icular set 



of abundan ces, chosen in our c ase to be the "coronal" set of 



fractions of 



Bryans et al. 



Feldmanl ( 119921 ). The ion 



( 120091 ) were used. This value and T e and emission measure (N%V) 
were set as free parameters to be determined, as well as those describing any nonthermal 
continuum present. A further component of the fitting function (drm_mod) allows for small 
modifications (gain and energy resolution) in the DRM. For this model fitting function 
(chianti_full), there is a total of 8 free parameters. There were slight disagreements 
in the fits, particularly around the energy of the Fe/Ni line complex (~ 8 keV), with the 
chianti spectrum underestimating the line flux, worsening the Xred- Figure Q] shows a 
RHESSI detector 4 Al spectrum and fit with the chianti_full function and normalized 
residuals. For this fit, 80 energy bins were used with 8 free parameters defining the model 
spectrum, giving Xmd = 1-07. With 72 degrees of freedom, this indicates a probability of 



Table 1 — Continued 



Flare Date RHESSI Time GOES CHIANTI_full Analysis CHiANTi_CONT Analysis 

Range (UT) Class Mean value A(Fe) A(Fc) 

Ar fu n oiF A(Fc) N cont Fe-Line Fe/Ni-Line 

Complex Complex 



11.. 


2003 


Apr 23 


01:00 - 


01:40 


M5.2 


43 


0.73 ±0.27 


7.96 ±0. 


14 


101 


7.99 ±0. 


.05 


8.07 ±0. 


.09 


12.. 


2003 


May 29 


01:10 - 


01:42 


XI. 1 




indeterminate 






118 


7.90 ±0 


.05 


7.99 ±0. 


.09 


13.. 


2003 


Aug 19 


10:00 - 


10:26 


M2.7 


121 


0.85 ±0.37 


8.03 ±0. 


16 


106 


7.94 ±0 


.05 


7.88 ±0. 


.12 


14.. 


2003 


Oct 22 


20:16 - 


20:37 


M9.9 




indeterminate 






23 


7.74 ±0 


.04 


7.83 ±0. 


.08 


15.. 


2003 


Oct 23 


20:06 - 


20:37 


XI. 1 


37 


0.58 ±0.14 


7.86 ±0, 


.10 


39 


7.94 ±0 


.04 


8.06 ±0. 


.09 


16.. 


2003 


Nov 2 


18:37 - 


18:59 


X8.3 


31 


0.50 ±0.06 


7.80 ±0. 


.05 


71 


7.88 ±0 


.08 


8.03 ±0. 


,12 


17.. 


2003 


Nov 11 


15:34 - 


16:33 


C8.5 


21 


0.64 ±0.17 


7.91 ±0. 


.10 


85 


8.00 ±0 


.04 


8.13 ±0. 


.09 


18.. 


2004 


Jan 5 


04:05 - 


04:52 


M6.9 


14 


0.53 ±0.05 


7.82 ±0. 


04 


37 


7.93 ±0 


.03 


8.08 ±0. 


.02 


19.. 


2004 


Jul 20 


12:40 - 


13:41 


M8.7 


41 


0.60 ±0.15 


7.88 ±0, 


10 


91 


7.92 ±0 


.06 


8.00 ±0. 


.10 


20.. 


2005 


Jan 16 


01:29 - 


02:27 


X2.6 


98 


0.28 ±0.06 


7.55 ±0, 


08 


102 


7.72 ±0 


.09 


7.90 ±0. 


.18 



a Af u n is the number of spectra analyzed with CHIANTI_full, N cont the number with CHIANTI.CONT. 

b "Indeterminate" : standard deviation in F > F or the standard deviation in A(Fe) > 0.3. See text for definition of F. 
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32% of exceeding this value through random statistical fluctuations in the count rates. The 
mismatc h at ~ 8 keV is unlikely to be due to the omission in chianti of n > 5 lines of 



Fe XXV ( [Phillips 



20081 ); more probably it is due to an instrumental line from the tungsten 
collimator grids not allowed for in the DRM. The fit gives an abundance factor F = 0.372, 
or A(Fe) = 7.67. 

A more general fitting function was also chosen that includes free-free and free-bound 
continua alone as calculated by chianti and line features with Gaussian profiles to fit 
the Fe and Fe/Ni line complexes. This fitting function (chianti_cont) has a single 
temperature and emission measure describing the continuum, and the two line fluxes as 
free parameters. The line energies were kept fixed at their theoretical values (6.65 keV and 
8 keV); there is only a slight temperature dependence of the line energies ~ 0.05 keV). As 
free-bound emission contributes to the continuu m, there i s some dependence on element 



Feldmanl ( 119921 ) abundance set for this 



abundances; as an initial estimate, we chose the 
calculation. As with the chianti_full function, gain and energy resolution adjustments 
with drm_mod and a nonthermal continuum were included. This gave a total of 9 free 
parameters. The observed and fitted spectra for the same interval and detector are shown 
in Figure El with normalized residuals. The value of Xred = 0-^0 implies a 89% probability 
of exceeding this value; there is a fairly random distribution of residuals with energy over 
the fitted range (5.67-33 keV). Allowing the 8-keV line flux to be a free parameter removes 
the enhanced residuals at this energy. 

As the Fe abundance will be derived from the line fluxes and the continuum values 
of T e and N^V, the dependence on element abundances of the continuum flux in the 
neighborhood of the Fe and Fe/Ni line complexes should be investigated, in particular the 



Fe abundance. Anticipating the discussion in Section El an Fe ab undance o 



or Fe/H = 8.13 x 10 is derived, i.e. a factor 1.55 less than the 



A(F e) = 7.91 



Feldmanl (119921 ) coronal 
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SPEX HESSI Count Rate vs Energy with Fit Function, Interval 23 
1 00.000 F r~i i ' i | i i 1 = 




-2 I. 



Energy (keV) 

Fig. 1. — (Upper panel) RHESSI detector 4 Al background-subtracted count-rate spectrum 
(black histogram) for an interval during the decay of the 2003 October 23 flare (energy range 
is 5 to 40 keV). The least-squares fitted spectrum (red) consists of two components folded 
through the DRM: (a) a thermal spectrum (vth, in green) including all lines and continua in 
CHIANTI; (b) a nonthermal photon spectrum (thick2, dark purple). The function (drm_mod) 
allows the energy resolution and gain to vary to give the best fit to the Fe line complex. The 
energy range of the fit was 5.7-33 keV (vertical dotted and dashed lines). The reduced 
chi-squared Xred = 1-07. The pink histogram is the background spectrum. (Lower panel) 
Residuals normalized to the statistical la uncertainties in each energy bin. 
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value. The effect on the total continuum of this abundance difference can be partly tested 
with RHESSI software, since a provision is made in the analysis software to adjust the 
abundances of Si, S, Ca, Fe, and Ni by different factors in the continuum function. We 
took spectra during the 2003 October 23 flare with the chianti_cont model function but 
reducing all the Si, S, Ca, Fe, and Ni abundances by a factor 1.55. The abundances for thes e 



Fludra fc Schmelz 



l ow-F IP elements would then be equal to the "hybrid" abundances by 
(119991 ). (There is nevertheless evidence that a constant reduction of all low-FIP elements 
i s not observed; estima tes of the flare potassium abundance from the RESIK instrument 



(jSylwester et al 



value, 5.67, of 



2010a) 



indica te A(K) = 5.86 ± 0.23 may be enhanced over the coronal 



Feldmanl (jl992j).) For the interval shown in Figure [2j the temperature was 



practically unchanged at 1.28 keV = 14.8 MK but the emission measure was 16% higher 
(0.57 x 10 49 cm -3 ). This difference folds directly into the derived Fe abundance. 

The contribution that Fe alone makes to the free-bound and total continuum is, 
however, much less than the total of Si, S, Ca, Fe, and Ni. The exact contribution was 
calculated at temperatures typical of those found in this analysis, the results being given in 
Table [2], where the percentage contribution of the Fe free-bound continuum at an energy of 
10 keV (near the peak of the RHESSI count rate spectrum in the Al state) is given to the 
total as a function of temperature (the contributions are similar at other energies). Thus, 
for T e = 25 MK, where the differences are most marked, the total continuum flux for a 
coronal Fe abundance (A(Fe) = 8.10) and an Fe abundance A(Fe) = 7.91 decreases by 6%, 
from 15 % to 21 %. This means that the derived emission me asure will incre ase by 6% if 



the other elements (Si, S, Ca, and Ni) remain at their coronal (lFeldmanlll992l ) abundances. 



At present, there is no recent detailed abunda nce analysis for these elements to confirm 



whether their abundances are the same as the 



Feldmanl (119921 ) abundances. However, it is 



worth pointing out that RESIK measureme nts of the continuum flux at somewhat lower 



energies (2.9-3.9 keV) (IPhillips et al 



20101 ) are better described by coronal abundances 



13 



( IFeldmanl Il992f ) than other abundance sets. Column 2 of Table [2] gives the percentage 
of the total free-bound continuum to the total continuum (free-free and free-bound) at 
10 keV. The calculations in Table [2] are confirmed by more detailed calculations involving 
all significant elements (J. Sylwester, work in preparation). 

Because of the mismatch using the chianti_full fitting function at ~ 8 keV, there were 
fewer spectral fits having small values of reduced Xred than with the chianti_cont fitting 
function. Temperatures and emission measures from chianti_cont and chianti_full 
are compared in the plot shown in Figure [3] for time intervals during three flares for which 
Xred < 1-5- It is evident that temperatures from the chianti_full model are smaller than 
those from chianti_cont by ~ 1 MK and emission measures larger by a factor ~ 2.5 (0.4 
in the logarithm). Although the specific reasons for these differences are unclear, they may 
be due to attempts in the fitting process with chianti_full to correct for the ~ 8 keV 
mismatch with a continuum function that is slightly too steep at energies > 8 keV. 

In our previous analysis, we compared spectral fits from seven of the nine RHESSI 
detectors suitable for spectral analysis in the low-energy region (detectors 2 and 7 have 
poorer energy resolution and higher energy threshold and so they were not used). Here, this 
comparison was done more systematically using a 30-minute time interval in the decline of 
the flare on 2002 July 26 with the chianti_cont fitting function. From a detailed analysis 
of the temperature and emission measure over this period, the mean la uncertainties (in 
MK) for the detectors 1, 3, 4, 5, 6, 8 and 9 were found to be respectively 1.20, 1.19, 1.05, 
3.02, 1.12, 1.08, 1.16. They are thus smallest for detector 4 and largest for detector 5. This 
is also true for the uncertainties in the emission measure estimates. In light of these results, 
our choice of spectral fits from detector 4 in our earlier work appears to be vindicated, and 
so we chose to fit spectra from this detector. 
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SPEX HESSI Count Rate vs Energy with Eit Function, Interval 23 
100.000 



I 0.000 
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0.010 
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d 



rval 23 Oh, 
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0.487,1 .27 
1 .70e + 004 
1 .62e+003 
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mod 0.85 



-I 1 14 



ick 
1.00 
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+ drm_mod 
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6.6^,0.100 
p. 100 

.6,150.,3.37,20.0,3.20e + 004 
,0.07 54,1.00 



n 



If 




Energy (keV) 



Fig. 2. — Same as Figure [T] with the thermal spectru m used for that plot 



repla ced with the 



20001 ) abundances 



thermal continuum given by CHIANTI with coronal (jFeldman fc Laming] 
and two lines with gaussian profiles at the mean energies of the Fe and Fe/Ni line complexes 
(6.65 and 8.0 keV respectively). The energy range is 5 to 40 keV. The value of Xred = 0.80. 



The mismatch in the residuals at ~ 8 keV has been significantly reduced. 



- 15 - 



20 



< 
x 
o 



1 6 



1 4 



12 



10 





10 



12 



14 



20 



1.0 



-0.6 



-0.2 



+ 0.2 



3 



.2 < 
x 
o 



.6 E 
o 

CD 

o 



1.0 



o 



+ 0.2 



Temp. (MK) from CHIANTLCONT log EM (10 49 cm" 3 ) from CHIANTLCONT 



Fig. 3. — Temperatures (MK) and logarithms of emission measures (unit = 10 49 cm -3 ) from 
spectral fits to 3 flares using fitting function chianti_full plotted against those from fitting 
function chianti_cont. The three flare intervals are numbers 2, 3, and 6a in Tabled], with 
symbols indicated in the legend. 
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2.4. Use and Validation of CHIANTI Spectra 



The RHESSI analysis code now uses version 6 of the chianti atomic code and 



database (IDere et al. 



20091 ) which includes the best atomic data available for the lines in 
the RHESSI range discussed here. In particular, an abundance of Fe is directly determined 
from comparison of the Fe line and Fe/Ni line complex fluxes with that from chianti 
using T K and emissi on measure from the continuum rather than with our earlier work 



(jPhillips et al. 



20061 ) in which the equivalent width of each line complex was determined 



from RHESSI spectra a nd comp a red w ith calculations based on the sum of individual 



lines in each complex by 



Phillips! (12004 ). Expressed as contribution functions Gp e (T e ) and 



Cpe/NiC^e)) or the amount of emission from the Fe line complex (defined to be all lines 
in the 6.391-7.005 keV range) and Fe/Ni line complex (all lines in the 7.728-8.907 keV 



range) from a plasma wit h unit vo 
few per cent between the 



ume emission measure, there are differences of only a 



Phillips! (12004 ) and chianti v. 6 calculations for the Fe line. 
There are more significant differences for the Fe/Ni line for T e < 12 MK, for which the 
chianti v. 6 calculations are higher by amounts that depend on T e . These are due to the 
addition of Fe XXIV satellites near the Fe XXV Is 2 — lsnp (n = 4, 5) lines not included by 
Phillips! (12004 ) or earlier versions of chianti. The chianti v. 6 calculations for GF e /m(T e ) 



are still incomplete in that Fe XXV Is 2 — lsn p (n > 5 ) lines and the associated satellites 



are not included. However, based on work by iPhillipsi (120081 ). this is unlikely to lead to 
an underestimate of Gpe/Ni(7e) by more than a few percent in the range of temperatures 
(approximately 10-22 MK) considered in this work. 

Validation of chianti spectra is possible by comparing with solar flare spectra from 
the Solar Maximum Mission Bent Crystal Spectrometer for the Fe line complex. There 
are small differences apparent in the ratio of some of the intense Fe XXIV satellites to the 
Fe XXV resonance line which provides the means of determining T e in high-resolution, 
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crystal spectrometer data. The BCS spectra analyzed by iLemen et all ( 119841 ) give 
T e = 15.0 ± 0.5 MK, whereas a re-analysis with chianti v. 6 leads to a higher temperature, 
16.5 MK. The difference is probably due to improved atomic data for both the Fe XXIV 
satellites and the collisional rates for the Fe XXV resonance line. Rather large differences 
(up to 0.01 keV) are apparent in the energies of Fe XXI-Fe XXIII satellites, though these 
are not likely to affect the total flux of the Fe line complex. There are no high-resolution 
flare spectra of the Fe/Ni line complex, so validation is not possible. 



Fe ABUNDANCE ANALYSIS 



Applying the chianti_full fitting function to RHESSI spectra gives temperature, 
emission measure, and the abundance factor F, which is the factor applied to the baseline 
Fe and Ni abundance and is determined by the fluxes of the Fe and Fe/Ni line com plexes. 



Feldman fe Lamingl ( 120001 ) (for Fe 



The baseline abundance set used was the coronal one of 
this is A(Fe)pL = 8.10). For all the flares analyzed, the values of F are practically constant 
with time and are unrelated to variations in T e . Since the Fe and Fe/Ni line features in 
the energy range chosen are nearly entirely due to iron , the measured Fe abundance is 



20001 ) value, or in logarithmic terms, 



F multiplied by the coronal (IFeldman fe Lamingl 
A(Fe) meas = log 10 F + A(Fe) FL . Table [1] gives the mean value of F and standard deviation 
for each of the time periods listed with the number of spectra used to derive F having 
Xred < 1-5- For a few flares there were too few good-quality spectra to give a reliable value 
of F (marked "indeterminate"). For the remaining 18 sets of spectra, the mean value of F 
is 0.55 ± 0.08 (s.d.), and so A(Fe) meas = 7.90 ± 0.02. This is several standard deviations less 
than A(Fe) FLl despite the fact that the value of F is poorly d etermined for some fl ares. 



It is a factor 2.5 more than the photospheric Fe abundance of 



Asplund et al 



f l2009h . To 



determine the abundance of Fe from the chianti_cont model, we followed the procedure 
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for analyzing spectra from the RESIK crystal spectrometer (jSylwester et al.ll2010al Jbl). For 
a spectrum in a particular (zth) RHESSI time interval, the iron abundance is determined 
from the flux 5i of either the Fe line or the Fe/Ni line complexes using 



fi(Fe) 



G(T i )EM l 



(1) 



where the temperature and emission measure EMi are output from the zth best-fit model 
fun ction. Thus, while like F in the chianti_full fitting model, is a factor multiplying 



the 



Feldman &: Laming 



(120001 ) Fe abundance to give the measured Fe abundance, it 
is defined in terms of the measured fluxes of either the Fe or Fe/Ni line complex, the 
continuum temperature Tj and emission measure EMi, and the contribution function 
GpeiTi) or G Fe / Ni (Ti). 

The estimated fluxes in the Fe line and Fe/Ni line complexes divided by the continuum 
emission measure when plotted against continuum temperature can be compared with 
the theoretical G-F e (T e ) and Grpe/Ni(^e) functions calculated with A(Fe)pL — 8.10 and the 
photospheric abundance (A(Fe) = 7.50). This is done in Figure H] (left) for the 2003 
August 19 flare. The uncertainties (from the OSPEX software) are larger for the Fe/Ni line 
complex, which for the temperatures analyzed here is about a factor 10-20 less intense than 
the Fe line complex, and are larger for lower-temperature spectra late in the flare when the 
emission was weaker. 

Values of measured Fe abundances from Eq. (1) were obtained for each spectral 
interval in this and other flares, and histogram distributions found for the Fe and Fe/Ni line 
complexes. This is illustrated for the Fe line complex in Figure H] (right) where numbers of 
values falling into intervals of 0.05 in A(Fe) are plotted, together with a best-fit gaussian 
curve. From the best-fit curves, a mean value of A(Fe) and standard deviation can be 
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Table 2. Percentage of Total Continuum at 10 keV due to Fe Free-bound 

Radiation 



Temperature (MK) 


Ratio free— bound to 
total continuum (%) a 


Photosphcric 
abundance (%) 


2.6x phot, abund. 

(%) 


Coronal abundance 
(4x photosph.) (%) 


10 


77 


4 


5 


8 


12 


73 


6 


8 


12 


15 


67 


8 


12 


18 


20 


58 


9 


15 


21 


25 


50 


8 


15 


21 



a Estimated with coronal I Feldman fe Lamindl2000r ) abundances (for which the Fe abundance is 4 times photo- 
sphcric.) 



A(Fe) from 2003 Aug. 19 (1000-1024) spectra 




18 20 22 '-° 

Temperature (MK) A(Fe) [A(H) = 12] 



Fig. 4. — (Left) Line flux/emission measure (10 49 cm" 3 ) for the Fe (6.65 keV) and Fe/Ni (8 
keV) line complexes vs. T e derived from RHESSI Al spectra during the decline of the flare 
on 2003 August 19. The GF e (T e ) and GrFe/Ni(^e) curves calculated from chianti are shown 
for coronal and photospheric Fe abundances. (Right) Distribution of the Fe abundances 
derived from the Fe line fluxes (histogram), with coronal and photospheric Fe abundances 
shown. The best-fit gaussian (smooth curve) indicates a mean value ^4(Fe) = 7.94 ± 0.05. 
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derived. The values derived are listed for each of the time periods in the last two columns 
of Table[TJ For some time periods, the Fe/Ni line complex was too weak to be measured 
so the Fe abundance could not be derived; the abundances are marked "indeterminate" 
when the standard deviation in A(Fe) was larger than 0.3 (corresponding to a factor 2 
in the abundance estimate). Similar plots to Figures H] were done for the time ranges in 
Table [TJ The total of all flares for the Fe line and Fe/Ni line complexes is shown in Figure 
There is a clear trend for the values of flux divided by emission measure to be below the 



coronal G(T e ) curves by a constant factor for the Fe line complex, s uggesting (as with t 



lie 



CHIANTI_full analysis) that the Fc abundance is smaller than the 



Feldman fe Laming 



(120001 ) value, but larger than the photospheric value. For the Fe/Ni line complex, the trend 
is similar at high temperatures (T e > 16.5 MK, but the points become steadily higher than 
the coronal abundance curve for decreasing temperatures. This may be because of the 
instrumental line at ~ 8 keV mentioned earlier. 

Figure [6] shows the distribution of Fe abundance values for spectra during all flares 
lumped together for the Fe line and Fe/Ni line complexes, with best-fit gaussian curves. 
The mean value of A(Fe) is 7.91 ± 0.10 from the Fe line, 8.01 ± 0.16 from the Fe/Ni line 
complex. The larger Fe abundance and larger uncertainty from the Fe/Ni line reflects the 
departure of the points from the theoretical GFe/Ni(7e) curve, so the Fe abundance from the 
Fe line complex is clearly the preferred value. We note that the continuum in this analysis 
has a small contribution from F e emission , so th ere is a slight dependence on the coronal 
abundances used, which was the iFeldmanl ( 119921 ) set. Use of an Fe abundance A(Fe) = 7.91 



instead of ^(Fe)^ = 8.10, if the elements Si, S, Ca, and Ni remain at their coronal values, 
should lead to emission measures that are ~ 6% larger. Eq. (1) indicates that the derived 
value of fi and therefore the Fe abundance using the chianti_cont emission model with 
coronal abundances for the continuum will be too large by 6% (0.03 in the logarithm). The 
abundance variation from flare to flare, in spite of the factor-of-100 range in the GOES 
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X-ray importance and RHESSI count rates, is very small and well within the standard 
deviations for each flare. 

4. Discussion and Conclusions 

Our analysis of RHESSI solar flare spectra gives estimates of the Fe abundance 
from two different emission models, chianti_full (continuum and lines as given by the 
CHIANTI atomic code) and chianti_cont (continuum given by chianti, Fe and Fe/Ni 
line complexes at 6.65 keV and 8 keV separately fitted). Both methods give Fe abundance 
values that are constant to within the estimated uncertainties and higher than those 
derived from the photosphere or meteorites. The chianti_full analysis was found to 
suffer from poor fits to the spectral region around the Fe/Ni line complex at 8 keV, leading 
to generally higher values of Xred- Using only spectral fits with Xred < we found from 
18 time intervals that the mean A(Fe) to be 7.87 ± 0.02, the uncertainty being in the 
mean values for each of the 18 intervals. The chianti_cont analysis has the advantage 
that the Fe and Fe/Ni line complexes can be fitted separately and A(Fe) can be derived 
for each. The possible presence of an instrumental line at ~ 8 keV adds to the observed 
Fe/Ni line complex and so the measured flux is an upper limit to the solar value. This 
probably explains the 23% difference between the values A(Fe) = 7.91 ± 0.10 for the Fe line 
during 25 time intervals and A(Fe) = 8.01 ± 0.16 for the Fe/Ni line during 21 time periods, 
though the Fe/Ni line complex is also a much weaker feature. When all the chianti_cont 
observations are lumped together (Figures [5] and [6]), it is clear, particularly from the Fe 
line, that the Fe abundance is constant to within estimated uncertainties from flare to flare. 
The distribution of all the estimates from the Fe line leads to an Fe abundance given by 
A(Fe) = 7.91 ± 0.10, which is our preferred value. 

Our value is higher by a factor 2.6 ± 0.6 than Fe abundance estimates from the 



10' 




Temperature (MK) 



Fig. 5. — Line flux/emission measure (10 49 cm~ 3 ) plotted against T e for all flares in this 
analysis for the Fe (6.65 keV) (upper panel) and Fe/Ni (8 keV) (lower panel) line com- 
plexes derived from RHESSI Al spectra. The calculated GF e (T e ) and G-p e /miT^) curves are 
shown as solid lines for coronal (IFeldman fe Laming||2000l ) and dashed lines for photospheric 



(lAsplund et al. 



ines for coronal (Fe ldman fe Laming 

~ 1 
2009) Fe abundances. 
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Phot. | Corona 

Mean A(Fe): 8.01 +/- 0.16 




7.4 7.6 7.8 8.0 8.2 
A(Fe) [A(H) = 12] 



7.4 7.6 7.8 8.0 8.2 
A(Fe) [A(H) = 12] 



Fig. 6. — Distribution of values of derived Fe abundance values, A(Fe)i, from fits to RHESSI 
spectra using chianti_cont for all flares in this analysis, compared with best-fit Gaussian 
curves. Left: from the Fe (6.65 keV) line complex; right: from the Fe/Ni (8 keV) line 
complex. The coronal and photospheric Fe abundances are shown as dashed vertical lines. 
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photosphere (e.g. A(Fe) = 7.50 ± 0.04: lAsplund et al 



than those from meteorites (7.45 ± 0.01: 



Lodders et al. 



(2009) 



those given for coronal plasmas, w hich range from 7.65 ( [Parkinson 



and by a factor 2.9 ± 0.7 



2009) ). However, it is lower than 



1977 



Fludra fe Schmelz 



1999) to 8.19 (White et al. 



Feldman fe Laming! (120001 ) and 



2000); it is a facto r 1.55 ± 0.5 lower than that given by 



Feldmanl (119921 ). viz. A(Fe) = 8.10, which is used as the 



coronal abundance set in chianti. The prel iminary results for fla res observed with the 



XRS instrument on Mercury MESSEN GER ( 



that is a factor 2.3 times photospheric (IGrevesse fe Sauval 



Dennis fe Starr 



2008) give an Fe abundance 



19981 ). i.e. A(Fe) = 7.9, which 



is consiste nt with our preferr e d val ue. Our value is within la of the Fe abundance of 



Table 2 of 



Fludra fe Schmelzl ( 119991 ). A(Fe) = 7.83, for their "hybrid" abundance model, 



taking account of the uncertainties in both our value and the hybrid model. Our value of 
7.91 ± 0.10 is in very close agreement (0.02 less) with the energetic particle abundances 



reported by iReamesI (119951 ) for gradual events, viz. A(Fe) = 7.93 ± 0.01, suggesting a 



relationship of the fractionation process involved in the formation of hot flare plasmas and 
the acceleration of solar energetic particles. The constancy of the RHESSI Fe abundance 
estimates points to fractionation processes in flares, at least their declining stages, to 
be similar to those f or the quiet Sun. It is rather against expectations, as discussed by 



Asplund et al. 



(120091 ). who state that "the degree of chemical separation varies significantly, 



being more severe in regions of higher solar activity." Though our measurements are during 
the decay of flares, the flare plasma is unlikely to be mixed with other plasma after the 
impulsive stage when (as is widely accepted) chromospheric evaporation occurs. 

If the Fe abundance derived here for flares is representative of the quiet solar corona 
and active-region levels, there are consequences for the radiation loss curve (Figure [7]). Iron 
ions are an important emitter for solar plasmas in the temperature range spanning values of 
the quiet corona to the tens of MK of solar flares. A maximum in radiation loss at ~ 1 MK 
is due to emission lines of Fe IX-Fe XV between 171 A and 284 A, and a second maximum at 
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~ 10 MK is due to Fe XVII-Fe XXIV X-ray lines. The precise abundance of iron is therefore 
important for studies of the energy balance in quiet coronal or flare loo ps, particularly 



for h igher densities when radiation loss may dominate conduction losses (IKlimchuk et al. 



20081 ). If the value obtained here, A(Fe) = 7.91, applies to all coronal plasmas with 
temperatures 1-20 MK, the radiation loss curve will be correspondingly modified. Figure [7] 
shows the radiation loss for optically thin plasma with this temperature range and for Fe 
abundances with photospheric and coronal values including our own. The contribution 
from Fe ions and hydrogen alone is also shown, with maxima at ~ 1 MK and ~ 10 MK. 
(The maximum at ~ 20000 K is due to Ly-a emission of hydrogen.) The radiation loss 
is less for a photospheric Fe abundance than a coronal Fe abundance by a factor ~ 3.4 
for quiet coronal loops (T e ~ 1 MK) and a factor ~ 3.0 for flare loops with T e ~ 10 MK. 
For A(Fe) = 7. 91 obtained in this work, the radiation loss is less than that for coronal 



Fe abundance (IFeldman 



19921 ) by a f actor 0.8 (X ~ 1 MK) and 0.72 (T e ~ 10 MK). The 



radiation loss curve approximated by 



Klimchuk et al 



(120081 ) with a piece -wise continu ous 



temperature- dependent form assumed abundances t hat are twice those of 



Meyer! (119851 ) gives A(Fe) = 7.6 for the corona, the 



Klimchuk et al. 



Meverl (119851 ): as 



(120081 ) value (7.9) is 



very nearly the one obtained here. 

The various models advance d up to the mid-1990s for explaining the FIP effect have 



been reviewed by iHenouxi (Il998l ). Fractionation according to the first ionization potential 
is generally explained by the fact that low-FIP (FIP < 10 eV) elements are partly ionized 
in the photosphere but high-FIP elements are predominantly neutral. It is difficult to 
obtain from any of these early models a quantitative enhancement of low-FIP elements 
in coronal plasmas, however, and it appears to be a feature of the present work that for 
flares with a range of X-ray importance that the enhancement of Fe is a particular value, 
constant from flare to flare. The constancy of the Fe abundance in flare plasmas is also 
consistent with the analysis of 2795 RESIK spectra indicating a constant enhancement (of 
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Radiative loss curves 




10 4 10 5 1 6 10 7 10 s 10 9 

Temperature (K) 



7. — Radiation los s curves for photospheric ( lAspfund et al.l 120091 ) and coronal 
( Feldman fc Lamingl 120001 ) abundances, and coronal abundances with the Fe abundance 
from this work (A(Fe) = 7.91). 
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a factor 5 over the photospheric abundance) of potassium in flares f jSylwester et al.ll2010af ) 
and an argon abundan ce that is also constant (in this case to within 20% of photospheric 



Ar abundance proxi es (jSylwester et al 



([Laming 



2004 . 



2010bl )). The more recent FIP model of Laming 



20091 ) involves the ponderomotive force that arises when Alfven waves pass 
from the chromosphere to the corona, and definite predictions about the enhancements of 
various elements can be made according to the dimensions of the coronal loop that the 
waves are incident on, its magnetic field, and the wave intensity. A particular example is 
given for a loop with length of 100,000 km and field strength of 7 G. RHESSI flare loops are 
likely to be much shorter, but it is interesting that for a fairly wide range of wave energy 
fluxes the coronal enhancement of Fe is between 2 and 3, as obtained in this work. 

In a wider context, it has been found that the occurrence of giant planets around sta rs 



depends sensitive 


y on 


Valenti & Fischer 


2008) 



y on metallicity, specifically the Fe abundance ( iFischer fe Valenti 



2005 



20081 ). This refers to Fe abundances in the stellar photospheres. However, 
it would be interesting to use the methods of this work to derive the coronal or flare 
abundance of Fe from the 6.65 keV line feature and compare with photospheric abundances 
to see if there are correlations of the difference with the probability of planet formation. 
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